Single-nucleon potential decomposition of the nuclear symmetry energy 
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The nuclear symmetry energy Esym.{p) and its density slope L{p) can be decomposed analyti- 
cally in terms of the single-nucleon potential in isospin asymmetric nuclear matter. Using three 
popular nuclear effective interaction models which have been extensively used in nuclear structure 
and reaction studies, namely, the isospin and momentum dependent MDI interaction model, the 
Skyrme Hartree-Fock approach and the Gogny Hartree-Fock approach, we analyze the contribution 
of different terms in the single-nucleon potential to the Esym{p) and L{p). Our results show that the 
observed different density behaviors of Esym{p) for different interactions are essentially due to the 
variation of the symmetry potential Uaym,i{p,k). Furthermore, we find that the contribution of the 
second-order symmetry potential Usym,2{p, k) to the L{p) generally cannot be neglected. Moreover, 
our results demonstrate that the magnitude of the Usym,2{p,k) is generally comparable with that 
of Usym,i{p,k), indicating the second-order symmetry potential Uaym,2{p,k) may have significant 
corrections to the widely used Lane approximation to the single-nucleon potential in extremely 
neutron(proton)-rich nuclear matter. 

PACS numbers: 21.65.Ef, 21.30.Fe, 21.60.Jz 
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I. INTRODUCTION 



During the last decade, the nuclear symmetry energy 
Esym{p) that essentially characterizes the isospin depen- 
dent part of the equation of state (EOS) of asymmetric 
nuclear matter has attracted much attention from dif- 
ferent fields due to its multifaceted influences in nuclear 
physics and astrophysics as well as some interesting 
issues regarding possible new physics beyond the stan- 
dard model [7Hl(j|. For example, the density slope L of 
the symmetry energy at nuclear matter saturation den- 
sity po has been shown to be important in determining 
several critical quantities such as the size of neutron- 
skin in heavy nuclei [lTl - [20j , location of neutron dripline 
dH, core-crust transition density 0, i, [H [H-il] and 
gravitational binding energy [2^ of neutron stars. The 
symmetry energy may also have significant influence on 
gravitational wave emission from compact stars [27l - l3l| . 
Furthermore, knowledge on the symmetry energy might 
be useful to understand the non-Newtonian gravity pro- 
posed in grand unified theories and to constrain prop- 
erties of the neutral weakly coupled light spin-1 gauge 
f7-boson originated from supersymmetric extensions of 
the standard model [l^, fs^ - lssj . 

In recent years, a great deal of experimental and theo- 
retical efforts have been devoted to determining the den- 
sity dependence of the symmetry energy (B], @ . Although 
significant progress has been made, large uncertainties on 
Esym{p) still exist even around nuclear matter saturation 
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density, e.g., while the value of Esym{po) is determined 
to be around 30 ± 4 MeV, mostly from analyzing nu- 
clear masses, the extracted density slope L scatters in 
a very large range from about 20 to 115 MeV, depend- 
ing on the observables and methods used in the studies 
13, m [13 (See, e.g., Refs. [H for a review of 

recent progress) . To reduce the uncertainties on the con- 
straints of Esy„i{po) and L is thus of critical importance 
and remains a big challenge in the community. 

So far, information on Esym{po) and its density slope L 
is essentially obtained from theoretical model analysis on 
the experimental data of heavy ion collisions 

[1111,111- 

[3 , nuclear mass [H [i^ , excitation energies of isobaric 
analog stat es |47| pygmy dipole resonance of neutron- 
rich nuclei |48l. [49l| . isovector giant dipole resonance of 
neutron-rich nuclei (50l . [Flj , and neutron skin thickness 
[13, [l3l- these theoretical models, an energy density 
functional with a number of parameters is usually as- 
sumed apriori, and the model parameters are then ob- 
tained from fitting experimental data and the empiri- 
cal values of some physical quantities. Information on 
Esymipo) and L is then extracted based on the obtained 
model parameters. Since all the phenomena/observables 
arc in someway at least indirectly and qualitatively re- 
lated to the Esymipo) and L, it is very useful to directly 
decompose Esym{po) and L in terms of some relevant 
parts of the commonly used underlying nuclear effective 
interaction [s^]. This decomposition of the Esym{po) and 
L provides an important and physically more transpar- 
ent approach to extract information about isospin depen- 
dence of strong interaction in nuclear medium from ex- 
periments and understand why the predicted symmetry 
energy from various models is so uncertain [53j . 
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In a recent work |3a|. based on the Hugenholtz-Van 
Hove (HVH) theorem |54l. [ssf , it was indeed shown that 
both Esymipo) and L are completely determined by the 
single-nucleon potential in asymmetric nuclear matter 
which can be extracted from the nucleon global optical 
model potentials. In that work, the Lane approximation 
[sgI ] to the single-nucleon potential in asymmetric nuclear 
matter, i.e., C/„/p(p, (5, fc) w Uoip.k) ± Usym,i{p,'k)5, has 
been assumed and also the momentum independent nu- 
cleon effective mass has been used, therefore, the higher- 
order effects such as the contributions from the second- 
order symmetry potential Usym,2{p,k) and the momen- 
tum dependence of the nucleon effective mass which can 
also contribute to L have been neglected. So far, to our 
best knowledge, there is not any empirical or experimen- 
tal information or even any theoretical predictions on the 
second-order symmetry potential Usym,2{p,k). It is thus 
interesting and important to estimate the Usym,2{p,k) 
with some well-established theoretical models. 

The main motivation of the present work is to eval- 
uate the Usym,2{p,k) and estimate its contribution to 
L based on several popular nuclear effective interaction 
models which have been extensively used in nuclear struc- 
ture and reaction studies. Our results indicate that, al- 
though the momentum dependence of the nucleon ef- 
fective mass might not be important, the second-order 
symmetry potential Usym.2{p,k) might have nonnegligi- 
ble contribution to the L. Furthermore, we find that 
the magnitude of the Usym,2ip,k) is generally compa- 
rable with that of Usym,iip,k), indicating the second- 
order symmetry potential Usym.2ip,k) may have signifi- 
cant corrections to the widely used Lane approximation 
Un/pip,S,k) Uo{p,k) ± UsymAp,k)5 for the single- 
nucleon potential in extremely neutron(proton)-rich nu- 
clear matter, e.g., in neutron stars and the neutron- 
skin region of heavy nuclei. These results imply that 
it is important to extract experimentally information on 

Usym,2{p, k). 

The paper is organized as follows. In Section |lTl we 
briefly recall the definition of the symmetry energy and 
the symmetry potential in asymmetric nuclear matter, 
and then derive the explicit expressions for the single- 
nucleon potential decomposition of the symmetry energy 
and its density slope. The results and discussions are pre- 
sented in Section lnil A summary is then given in Section 
IIVI For completeness, the theoretical models used in the 
present paper are briefly described in the Appendix. 



II. THEORETICAL FORMULISM 

A. The symmetry energy and the symmetry 
potential in asymmetric nuclear matter 

Due to the exchange symmetry between protons 
and neutrons in nuclear matter when one neglects the 
Coulomb interaction and assumes the charge symmetry 
of nuclear forces, the EOS of isospin asymmetric nuclear 



matter, defined by its binding energy per nucleon, can be 
expanded as a power series of even-order terms in isospin 
asymmetry 5 as 



E{p,5) = Eq{p) + E,yn,{p)5^ + 0{S^) 



(1) 



where p = Pn + Pp is the baryon density and S ~ {p„ — 
Pp)/p is the isospin asymmetry with p„ and pp denoting 
the neutron and proton densities, respectively; Eq{p) — 
E{p, (5 = 0) is the EOS of symmetric nuclear matter, and 
the nuclear symmetry energy is expressed as 



Esymip) 



1 d^E{p, S)^ 
2! W 



\s=o- 



(2) 



The higher-order terms of S in Eq. ([T|) are negligible, 
leading to the well-known empirical parabolic law for the 
EOS of asymmetric nuclear matter, which has been ver- 
ified by all many-body theories to date, at least for den- 
sities up to moderate values (See, e.g., Ref. Q). 

Around the nuclear matter saturation density po, the 
nuclear symmetry energy Esymip) can be expanded as 



Esymip) = £^sym(po) + + O(x'), 



(3) 



where X — iP ~ Po)/3po is a- dimensionless variable and 
L = L(po) is the density slope parameter of the symmetry 
energy at po- More generally, the slope parameter of the 
symmetry energy at arbitrary density p is defined as 



dp 



(4) 



The slope parameter L at po characterizes the density de- 
pendence of the nuclear symmetry energy around nuclear 
matter saturation density po, and thus carry important 
information on the properties of nuclear symmetry en- 
ergy at both high and low densities. 

The single-nucleon potential Ur{p,S,k) (we assume 
T = 1 for neutrons and —1 for protons in this work) 
in asymmetric nuclear matter generally depends on the 
baryon density p, the isospin asymmetry 5 and the am- 
plitude of the nucleon momentum k. Due to the isospin 
symmetry of nuclear interactions under the exchange 
of neutrons and protons, the single-nucleon potential 
Ur{p, S, k) can be expanded as a power series of 5 as [52| 

Ur{p,5,k) = c/o(p,fc)+ UsymAp,k){Tsy 

i=l,2,--- 

= Uo{p,k) +UsymAP^k){TS) 



+ Usy,n,2{P^k){T5f +■ 



(5) 



where Uo{p,k) = t/„(p, 0, fc) = [/p(p, 0,fc) is the single- 
nucleon potential in symmetric nuclear matter and 
Usym.iip, k) are expressed as 



l 9-[/„(p,J,fc) , 

UsymAP^k) = — \s=0 

i-iyd'Up{p, S,k). 



1 5=0, 



(6) 
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with Usym^ip, k) being the well-known nuclear symmetry 
potential [Q| (where Usym.i is denoted by Usym), and the 
higher-order term Usym,2{p, k) being called as the second- 
order nuclear symmetry potential here. Neglecting the 
higher-order terms {6"^, S^, ■ ■ ■) in Eq. (O leads to the 
well-known Lane potential [sgI ]. i.e., 

Urip, S, k) « [/o(p, k) + UsyrnAP, k^rS), (7) 

which has been extensively used to approximate the 
single- nucleon potential Urip,S^k) in asymmetric nu- 
clear matter, and in this case the symmetry potential 

U sym 

k) can be obtained approximately by |6|,|57| 



U. 



sym, 1 



Un{p,S,k) - Up{p, S, k) 



26 



(8) 



B. Single-nucleon potential decomposition of the 
symmetry energy and its density slope 

According to the HVH theorem [s^, , the chemical 
potentials of neutrons and protons in asymmetric nuclear 
matter with energy density e{p, S) = pE[p, 5) can be ex- 
pressed, respectively, as 

t{kFj + UAp,5,kpJ = (9) 



tikF^) + Upip,S, kpj 



dpn 

de{p,S) 
dpp 



(10) 



where t(kp^) = kp /2m is the nucleon kinetic energy 

at Fermi momentum kp^ = kpil + tS^^^ with kp = 
(377^^/2)^/"^ being the Fermi momentum in symmetric 
nuclear matter at density p. We would like to point out 
that the HVH theorem is independent of the detailed na- 
ture of the nucleon interactions used and has been strictly 
proven to be valid for any interacting self-bound infinite 
Fermi system [s^, ■ 

The right-hand side of Eq. © can be further written 

as 



de{p,S) de{p,S) dp de{p,S) 85 
dpn 



dp dpn 06 dpn 
de{p,5) lde{p,5) de{p,6) 5 



dp 



p 85 



dS p 



(11) 



Similarly, the right-hand side of Eq. (|T0|) can be ex- 
pressed as 



de[p,5) 8e{p,6) lde{p,S) de{p,S) S 



dpp 



dp p dS 



d5 p 



(12) 



Subtracting Eq. from Eq. ([TT|) and noting e{p, 6) = 
pE[p, 5), we then obtain 



de{p,S) de{p,S) 



dpn 



dpp 



2de{p,S) _^dE{p,S) 



p dS 



d6 



(13) 



while adding Eq. (fTTj) and Eq. (|T2)) . we have 
de{p,S) de{p,5) 



dpn 
2E{p,6) + 2p 



dpp 

dE{p, S) 
dp 
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dE{p,6) 
dS ■ 



(14) 



On the one hand, substituting Eq. ([T|) into Eq. (fT3)) and 
Eq. (ITil) respectively leads to following expressions, i.e., 

tikpj - tikp^) + Unip,S, kpj - Up{p, (5, kp^) 
= AEsyUp)5 + 0{5^) (15) 



and 



t{kpj + t{kp^) + C/„(p, 6, kpj + Up{p, 6, kp^) 
dEoip) 



2Eoip) + 2p- 



dp 



+ [-L{p) - 2E,y^{p) 
+ 0{5^). 



(16) 



On the other hand, tikp^) and UT{p,S^kp^) can be ex- 
panded as a power series of S, respectively, as 

t{kp^)=tikp) 



dt{k) 
+ -dl^^' 



hpirS) 



and 



fc2 dH{k) 2kp dt{k) 
T dk^ '''^ ~ ~9 aA^' 
+ 0{S'), 

Ur (p, 5, kp^ ) = Uo {p, kp) 
kp dUn{p,k) 
3 dk 

kp dUgyrr,,l{p, k) 

3 dk 

kl d^Uo{p,k) 



(17) 



Ifep + Usym,lip, kp) (tS) 



u. 



sym, 2 



{p, kp) 



1 

21 9 
0(<53). 



dk^ 



_ 2kp dUojp, k) 
9 dk '''^ 



(18) 



Substituting Eqs. ([17]) and (fT8|) into the left-hand sides 
of Eqs. (|15p and ([TS]) . and comparing the coefficients of 
the Ist-order 6 terms in both left and right hand sides, 
we then obtain 

p . ^ 1,, f , idm+uo{pM, , 

t^symyP) = -j^UsymSyP, Kp) + \kp ■ kp, 

(19) 

while comparing the coefficients of 2nd-order 5 terms in 
both sides leads to the following expression 
3 

L{P) = -jUsy,ns{P^kp) +Wsym,2{Pikp) 

, dUsy^.i, , , 1 9[t(fc) + £/o(p,(5)] , 

\kF ■ Kp + -; \kF ■ Kp 



dk ' 6 

i 9^[t(fc) + [/o(p,j); 

6 dk"^ 



dk 



kp Kp. 



(20) 
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It should be stressed that higher-order terms of the 
smgle-nucleon potential in Eq. ([5]) (i.e., J7sym,3(p, fc) and 
higher-order terms) have no contributions to Esym{p) 



and and thus Eq. p9)) and Eq. pO|) are com- 

plete and exact, and are valid for arbitrary density p. 
Furthermore, Eq. (fT9|) and Eq. (|20| can be rewritten as 



L{p) 



3 2to;5 

2 , 



3 2toS 



-U. 



syrn, 1 



(P, kp) 



dU. 



sym,l 



dk 



\kp ■ kp + 3Usym,2{p, kp), 



(21) 
(22) 



in terms of the nucleon effective mass mg(/5, fc) in sym- 
metric nuclear matter which is generally dependent of 
the density p and the nucleon momentum fc, i.e.. 



ml{p,k) 



m 



1 I m dUo{p,k) 
^ Wk dk 



due to the following relations 



d[t{k) + Uo{p, k)] 

dk 

d^[t{k)^Uo{p,5)] 



h^k. 



h?k dml 



For convenience, we can reexprcss Esym(p) and L(p), 
respectively, as 

Esyrnip) = Ei{p)+E2{p), (23) 

L{p) = Li{p) + L2ip)+L3ip)+Li{p) + L5ip), 



with 



Eiip) 



1 h^kj. 



3 2m^{p, kp) 



E2{P) = -^Usym.liP'kp) 



Li{p) 



L2{P) 



2 



3 27775(/9, kp) 
1 ti'kl 



dml{p,k) 



Qml^{p,kp) dk 



\kF 



Ez{p) = -^Usymsip, kp) 



Li{p) 
L^ip) 



dU. 



sym,l 



iP,k) 



dk 



Ikp ■ kf 



3U. 



sy7n,2 



{P,kp). 



(24) 
(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



In this way, one can see that Ei (p) represents the kinetic 
energy part (including the effective mass contribution) of 
the symmetry energy while E2 (p) is due to the symmetry 
potential contribution to the symmetry energy. Further- 
more, Li{p), L2{p), L^ip), L4{p), and L^ip) have re- 
spective physical meaning, namely, Li{p) represents the 



kinetic energy part (including the effective mass contri- 
bution) of the L parameter, L2{p) is from the momen- 
tum dependence of the nucleon effective mass, Lz{p) is 
due to the symmetry potential contribution, L^i^p) comes 
from the momentum dependence of the symmetry po- 
tential, and L^^p) is from the second-order symmetry 
potential Usym,2ip, kp). In this way, the symmetry en- 
ergy Esym{p) and its slope L{p) have been decomposed in 
terms of Uo{p,k), Usyrn,i{p, k), Usym^P^k) and/or their 
1st and 2nd-ordcr partial derivatives with respect to k 
and 5. In particular, at nuclear matter saturation density 
Po, UoiPOik), Usym.i{po,k) and Usym.APo.k) (and thus 
Egymipo) and its slope parameter L) can be determined 
completely from the isospin and momentum dependent 
nucleon global optical potential, which can be directly ex- 
tracted from nucleon-nucleus and (p,n) charge-exchange 
reactions (See, e.g., Refs. [H, [H, [53). 

For the single-nuclcon potential decomposition of the 
slope parameter L{p), if we use the Lane potential Eq. 
([7]) and neglect the contributions from the momentum 
dependence of the nucleon effective mass, the L(p) is then 
reduced to 



L{p) ^Li + L3 + L4 



(31) 



namely, 



Lip) 



27,2 



3 2777?; 



+ -jUsyniAP^ kp) 



+ 



dUsymAp, k) 

dk 



Uf • ki 



(32) 



which has been used in previous work [38f . Although the 
Lane potential could be a good approximation in evaluat- 
ing Usym,i{p, k) as in Eq. ([8]) @,[53l, it would be interest- 
ing to see if the higher-order Usym,2{p,kF) contribution 
to the L{p) is significant or not. Using three popular nu- 
clear effective interaction models, we will demonstrate 
in the following that the contribution of higher-order 
Usym.2{p, k) term to L{p) generally cannot be neglected. 



III. RESULTS AND DISCUSSIONS 

In the following, we analyze the single-nucleon poten- 
tial decomposition of Esym{p) and L{p) as well as the 
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density and momentum dependence of Usym.iip,k) and 
Usym,2{p,k) using three popular nuclear effective inter- 
action models which have been extensively used in nu- 
clear structure and reaction studies, namely, the isospin 
and momentum dependent MDI interaction model, the 
Skyrme Hartree-Fock approach, and the Gogny Hartree- 
Fock approach. One can find the details about these 
three models in the Appendix. A very useful feature of 
these models is that analytical expressions for many in- 
teresting physical quantities, such as the single-nucleon 
potential in asymmetric nuclear matter at zero temper- 
ature, can be obtained, and this makes our analysis and 
calculations physically transparent and very convenient. 

For the MDI interaction, we use 3 parameter sets, i.e., 
X = —1, X = and x = 1 [36[, which give three different 
density dependence of the symmetry energy, namely, stiff 
one, moderate one, and soft one, respectively, and have 
been applied extensively in transport model simulations 
for heavy ion collisions. For the Skyrme interaction, we 
mainly use the famous SKM* l&ol a nd SLy4 [6l| as well 
as the recently developed MSLO [1^ . In addition, a num- 
ber of other Skyrme interactions are used for the single- 
nucleon potential decomposition of Esym{p) and L{p) at 
Pf). For the Gogny interaction, we use the existing Dl 
[ia], DIS [si, DIN m, and DIM [HH which have been 
successfully applied in nuclear structure studies. 



A. Single-nucleon potential decomposition of 
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FIG. 1: 

interaction with x - 



(Color online) Density dependence of Esym{p), 
and E2(p) = iUsy,n,i{p,kF) in the MDI 
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FIG. 2: (Color online) Same as Fig. [T] but in the Skyrme 
Hartree-Fock approach with MSLO (a), SLy4 (b), and SKM* 
(c). 
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FIG. 3: (Color online) Same as Fig. [T]but in Gogny Hartree- 
Fock approach with Dl (a), DIS (b), DIN (c), and DIM (d). 



similar behaviors can also be seen from Fig. [2] for the 
Skyrme-Hartree-Fock calculations and Fig. [3] for the 
Gogny-Hartree-Fock calculations. 

Furthermore, one can see from Fig. [TJ Fig. [5] and 
Fig. [3] that for all the interactions considered here, the 
Ei{p) increases with the density and is always positive 
while E2{p) can increase or decrease with the density 
and even become negative at higher densities. These fea- 



In Fig. [U Fig. [2]and Fig. [3l we plot the density depen- 
dence of Esymip), Ei{p) and E2{p) in the MDI interac- tures can be understood since the Ei{p) 



tion model, the Skyrme-Hartree-Fock approach, and the 
Gogny-Hartree-Fock approach, respectively. One can see 
from Fig. [l]tliat for different x values, the Ei{p) displays 
the same density dependence while the i?2(p) exhibits 
very different density behaviors, indicating that the dif- 
ferent density dependences oi Egymip) for a; = — 1, and 
1 are completely due to the different density dependence 
of i?2(p), i.e., the symmetry potential Ugym.iip, kp). The 



3 2raQ{p,kF) 

is determined uniquely by the single-nucleon potential 
Uo{p,k) in symmetric nuclear matter for which reliable 
information about its density and momentum depen- 
dence has already been obtained from heavy-ion col- 
lisions, see, e.g., Ref. 0, albeit there is still some 
room for further improvements, particularly at high mo- 
menta/densities, and the nuclear effective interactions 
arc usually constructed to describe reasonably Uo{p,k), 
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especially around pQ. However, on the contrary, the sym- 
metry potential Usym,iip, kp), which mainly reflects the 
isospin dependence of the nuclear effective interaction in 
nuclear medium, is still not very well determined, espe- 
cially at high densities and momenta. In fact, it has been 
identified as the key quantity responsible for the uncer- 
tain high density behavior of the symmetry energy as 
stressed in Ref. (See also Refs. ^,^)- These re- 
sults show that the observed different density behaviors 
of Esymip) for different interactions are essentially due 
to the variation of the symmetry potential J7sym,i(Pi k). 
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B. Single-nucleon potential decomposition of L(p) 

In order to illustrate the single-nucleon potential de- 
composition of we show in Fig. HI Fig. [5] and 
Fig. [5] the density dependence of L{p), Li{p), L2{p), 
L^{p)^ Li{p) and Lc,{p) in the MDI interaction model, the 
Skyrme-Hartree-Fock approach, and the Gogny-Hartree- 
Fock approach, respectively. It is seen that the Li{p) 
displays the similar density dependence for all the in- 
teractions considered here, just like the Ei{p) shown in 
Fig. m Fig. [Hand Fig. [3] because of Li{p) = 2Ei{p). 
The L2{p) is seen to contribute a small negative value 
to the L{p), indicating that the momentum dependence 
of the nuclcon effective mass is generally unimportant. 
In particular, one can sec from Fig. [5] that L2{p) van- 
ishes for the Skyrme-Hartrce-Fock calculations due to 
the fact that the nucleon effective mass is momentum 
independent for the zero-range Skyrme interaction con- 
sidered here. The L3(/o) exhibits different density depen- 
dences for different interactions, reflecting the variation 
of Usym,i{Pikp) with density for different interactions. 
The Li[p) represents the contribution of the momentum 
dependence of the symmetry potential to the L[p), and 
it displays different density dependence for different in- 
teractions and can be negative or positive, depending on 
the interaction used. 



FIG. 5: (Color online) Same as Fig. [4] but in the Skyrme 
Hartree-Fock approach with MSLO (a), SLy4 (b), and SKM* 
(c). 
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FIG. 6: (Color online) Same as Fig. |4]but in Gogny Hartree- 
Fock approach with Dl (a), DIS (b), DIN (c), and DIM (d). 
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FIG. 4: (Color online) Density dependence of I/(p), Li(p) = 



3 2mQ (p,kp) ' 
3^ 



L2{p) = 



6 ml^(p,kp) 



dk 



iKj^ - kp, and Ls{p) = 

sym,2{p, kp) in the MDI interaction with x = —1 (a), (b), 
and 1 (c). 



It is particularly interesting to analyze the L5 (p) since 
it reflects the higher-order Usym.2{p, kp) contribution to 
the L{p) and has been neglected in the previous work 
[sst . From Fig. HI Fig. [5] and Fig. [HI it is surprising 
to see that the L^[p) may play an important role to de- 
termine the L{p). In the MDI interaction with a; = — 1, 
the L^{p) is always positive and increases rapidly with 
density while the opposite behavior is observed for the 
MDI interaction with x = 1. For the MDI interaction 
with X ~ 0, the L^i^p) is positive and moderately in- 
creases with density. For the Skyrme-Hartree-Fock cal- 
culations, the Lq(p) can be positive or negative while it 
is always negative for the Gogny-Hartree-Fock calcula- 
tions with Dl, DIS, DIN, and DIM. These results in- 
dicate that generally the higher-order Usym.iip, kp) con- 
tribution to the L(p) cannot be neglected and the Lane 
potential approximation to the single-nucleon potential 
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in asymmetric nuclear matter may cause significant er- 
ror for the determination of L{p) from the single- nucleon 
potential decomposition. 

Due to the special interest on the values of Esym{p) 
and L{p) at poj we list in Table |T] the values of the char- 
acteristic parameters po: Eo{pq), Esym(po), Ei, E2, L, 
Li, L2, is, and L5 at po for the MDI interaction with 
X = —1, and 1, the Gogny-Hartrec-Fock predictions 
with Dl, DIS, DIN, and DIM as weU as the Skyrme- 
Hartree-Fock predictions with 112 standard Skyrme in- 
teractions (In the table, the interactions in different mod- 
els are in order according to the L value and the corre- 
sponding reference for different interactions is included 
in the last column). For the MDI interaction model, one 
can see from Table U that the value of E2{po) is com- 
parable with that of £'i(po), and all the Eo{po), Ei{po), 
and £'2(^0 ) are the same for different x values due to 
the fact that the Uo{p, k) and UsymA{po, k) are indepen- 
dent of the X parameter by construction [36| . Because 
of the same reason, the values of ii, L2, ^3, L4 are all 
independent of the x parameter too. Therefore, for the 
MDI interaction, the x dependence of the L parameter 
is completely determined by the second-order symmetry 
potential Usym,2ip, kp)- Depending on the value of the x 
parameter, the Usym,2{p, kp) can be positive or negative. 
In particular, we have L = Li + L2 + L3 + £4 = 48.5 
MeV if we assume Usym.2{p-, kp) = 0. Furthermore, it is 
seen that the L2 contribution is relatively small compared 
with that of Li, L3, or ^4, indicating that the contribu- 
tion due to the momentum dependence of the nucleon 
effective mass is unimportant, and this is consistent with 
the observation from Fig. 21 

For the Gogny interaction, the value of L listed in Ta- 
ble U ranges from about 18 MeV to 34 MeV. And simi- 
larly with the MDI interaction model, the L2 is relatively 
small, ranging from about —8 MeV to —3 MeV. It is in- 
teresting to see that the values of ii and L3 from the 
Gogny interactions are quite similar with those of the 
MDI interaction model, and thus the difference of the L 
parameter from different interactions in these two models 
is mainly due to the variation of the L4 and L5. Further- 
more, for different Gogny interactions considered here, 
the value of L4 can be positive or negative while the 
value of is always negative, and the L5 contribution 
to the L parameter usually is relatively important. 

For the standard Skyrme interactions, we have L2 = 
McV. For the 112 Skyrme interactions considered in 
Table HI it is seen that the value of L ranges from about 
-50 MeV to 160 MeV, Li from about 20 MeV to 60 MeV, 
L3 from about 4 MeV to 74 MeV, L4 from about -48 
MeV to 36 MeV, and L5 from about -102 MeV to 58 
MeV. Therefore, the different term contributions to the 
L parameter can change a lot in the standard Skyrme 
interactions, especially for L3, L4, and L5. 

In order to see more clearly and intuitively the differ- 
ent term contributions to the L parameter, we show in 
Fig. [7] the correlations of ii(po), -^3(^0), -^4(/Oo), and 
Lsipn) with L{po) for the MDI interaction with x — —1, 
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FIG. 7: (Color online) Correlations of Li{po) (a), Lz[po) (b), 
Liipo) (c), and L^i^po) (d) with L{po) for the MDI interaction 
with a; = — 1, and 1, Cogny-Hartree-Fock predictions with 
Dl, DIS, DIN and DIM as well as Skyrme-Hartree-Fock pre- 
dictions with the 112 standard Skyrme interactions considered 
in Table [J 



and 1, the Gogny- Hartree-Fock predictions with Dl, 
DIS, DIN, and DIM as well as the Skyrme-Hartree- 
Fock predictions with the 112 standard Skyrme inter- 
actions. One can see from Fig. [7] that the results from 
the MDI interaction model and the Gogny-Hartree-Fock 
calculations are essentially consistent with the system- 
atics of the Skyrme-Hartree-Fock predictions. Based on 
these calculated results, we find that a statistical analysis 
can lead to Li(po) « 30 ± 6.5 MeV, Lsipa) « 46 ± 9.5 
MeV, Liipo) « -4 ± 15 MeV, and £5(^0) ~ -35 ± 30 
MeV. These results indicate that, within the standard 
Skyrme-Hartree-Fock energy density functional, ii(po) 
and L^{po) are relatively well constrained, and the main 
uncertainties are due to the Li^po) and L^^po) contribu- 
tions. Furthermore, it is interesting to see from Fig. [7] 
that there exhibits an approximately linear correlation 
between L^^po) and L{pq). If we use the present empir- 
ical constraint L{po) = 60 ± 30 MeV, then we find the 
E^ipo) can vary from about —66 MeV to 24 MeV, i.e., the 
value of Usym.2{poy kp) can vary from about —22 MeV to 
8 MeV. 



C. The symmetry potential [/; 



si/m, 1 



{P,k) 



Shown in Fig. [5] is the momentum dependence of 
the Usym,i{p,k) at p = O.bpo, po and 2pQ using the 
MDI interaction with x = —1, 0, and 1. For compar- 
ison, we also include in Fig. [5] the corresponding re- 
sults from several microscopic approaches, including the 
rclativistic impulse approximation (RIA) [l03l Il04{ us- 
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TABLE I: The characteristic parameters po (fm"''), Eo{po) (MeV), Esy„^{po) (MeV), Ei (MeV), E2 (MeV), L (MeV), Li 
(MeV), L2 (MeV), L3 (MeV), L4 (MeV), and L5 (MeV) at saturation density po for the MDI interaction with a; = -1, and 
1, the Gogny-Hartree-Fock predictions with Dl, DIS, DIN and DIM as well as the Skyrme-Hartree-Fock predictions with 112 
standard Skyrme interactions. The interactions in different models are in order according to the L value. The corresponding 
reference is included as the last column. 



Model 


po 


Eoipo) 


Esym{po) 


El 


E2 


L 


Li 


L2 


Ls 


Li 


is 


Ref. 


MDI 


























MDI x=l 


0.160 


-16.2 


30.5 


18.2 


12.3 


14.6 


36.4 


-6.2 


36.9 


-18.6 


-33.9 


[36] 


MDI x=0 


0.160 


-16.2 


30.5 


18.2 


12.3 


60.2 


36.4 


-6.2 


36.9 


-18.6 


11.6 


[36] 


MDI x=-l 


0.160 


-16.2 


30.5 


18.2 


12.3 


105.7 


36.4 


-6.2 


36.9 


-18.6 


57.2 


[36] 


Gogny 


























Dl 


0.166 


-16.3 


30.7 


18.8 


11.9 


18.4 


37.7 


-5.0 


35.6 


-17.1 


-32.7 


[62] 


DIS 


0.163 


-16.0 


31.1 


17.9 


13.3 


22.4 


35.7 


-7.6 


39.8 


-7.5 


-37.9 


[63] 


DIM 


0.165 


-16.0 


28.6 


16.8 


11.8 


24.8 


33.6 


-4.2 


35.3 


4.5 


-44.3 


[65] 


DIN 


0.161 


-16.0 


29.6 


16.5 


13.1 


33.6 


33.0 


-3.7 


39.2 


9.2 


-44.2 


M 


Skyrme 


























Z-fit 


0.159 


-16.0 


26.8 


14.5 


12.3 


-49.7 


29.0 





36.9 


-13.3 


-102.4 


[67] 


Esigma-fit 


0.163 


-16.0 


26.4 


14.8 


11.6 


-36.9 


29.6 





34.9 


-18.1 


-83.2 




E-fit 


0.159 


-16.1 


27.7 


14.1 


13.6 


-31.3 


28.2 





40.7 


-15.2 


-84.9 


[67] 


Zsigma-fit 


0.163 


-15.9 


26.7 


15.9 


10.8 


-29.4 


31.8 





32.4 


-17.1 


-76.5 


[67] 


SVII 


0.143 


-15.8 


27.0 


11.4 


15.6 


-10.2 


22.8 





46.7 


-11.7 


-68.0 


[68] 


SkSC4o 


0.161 


-15.9 


27.0 


12.3 


14.7 


-9.7 


24.6 





44.0 





-78.3 


[69] 


SVI 


0.143 


-15.8 


26.9 


12.0 


14.9 


-7.3 


24.1 





44.6 


-12.2 


-63.8 


[68] 


ZsigmaS-fit 


0.162 


-16.0 


28.8 


16.0 


12.8 


-4.5 


32.0 





38.4 


-17.9 


-57.0 




v070 


0.158 


-15.8 


28.0 


11.6 


16.4 


-3.5 


23.1 





49.3 


-34.7 


-41.2 


[70] 


v075 


0.158 


-15.8 


28.0 


11.6 


16.4 


-0.3 


23.1 





49.3 


-27.8 


-45.0 


[70] 


SI 


0.155 


-16.0 


29.2 


13.2 


16.0 


1.2 


26.4 





48.1 


-11.4 


-61.9 


[67] 


v080 


0.157 


-15.8 


28.0 


11.6 


16.4 


2.2 


23.1 





49.3 


-21.7 


-48.5 


[70] 


v090 


0.157 


-15.8 


28.0 


11.6 


16.4 


5.0 


23.1 





49.3 


-11.6 


-55.8 


[70] 


Skz4 


0.160 


-16.0 


32.0 


17.5 


14.5 


5.8 


35.1 





43.4 


27.9 


-100.6 


\m 


SkSClS 


0.161 


-15.9 


28.0 


12.3 


15.7 


6.7 


24.6 





47.0 





-65.0 


[69] 


BSk3 


0.157 


-15.8 


27.9 


10.8 


17.1 


6.8 


21.6 





51.3 


-16.6 


-49.6 


[72] 


MSk3 


0.158 


-15.8 


28.0 


12.2 


15.8 


7.0 


24.3 





47.5 





-64.8 


[73] 


vl05 


0.157 


-15.8 


28.0 


11.6 


16.4 


7.1 


23.1 





49.3 





-65.3 


[70] 


MSk4 


0.157 


-15.8 


28.0 


11.6 


16.4 


7.2 


23.1 





49.3 





-65.2 


[73] 


BSkl 


0.157 


-15.8 


27.8 


11.6 


16.2 


7.2 


23.1 





48.7 





-64.7 


m 


vllO 


0.157 


-15.8 


28.0 


11.6 


16.4 


7.5 


23.1 





49.3 


3.2 


-68.1 


[70] 


MSk5 


0.157 


-15.8 


28.0 


11.6 


16.4 


7.6 


23.1 





49.3 





-64.9 


[73] 


BSk2p 


0.157 


-15.8 


28.0 


11.6 


16.4 


7.8 


23.2 





49.3 


-14.9 


-49.7 


[75] 


BSk2 


0.157 


-15.8 


28.0 


11.7 


16.3 


8.0 


23.3 





49.0 


-14.8 


-49.6 


[75] 


MSk8 


0.158 


-15.8 


27.9 


11.0 


16.9 


8.3 


22.1 





50.6 





-64.5 


[76] 


vlOO 


0.157 


-15.8 


28.0 


11.6 


16.4 


8.7 


23.1 





49.3 


-3.5 


-60.2 


[70] 


MSk7 


0.158 


-15.8 


27.9 


11.6 


16.4 


9.4 


23.1 





49.1 





-62.9 


[77j 


MSk6 


0.157 


-15.8 


28.0 


11.6 


16.4 


9.6 


23.1 





49.3 





-62.8 


[73] 


SIII 


0.145 


-15.9 


28.2 


15.1 


13.1 


9.9 


30.2 





39.2 


-14.8 


-44.7 


[78] 


MSk9 


0.158 


-15.8 


28.0 


12.2 


15.8 


10.4 


24.3 





47.5 


0.0 


-61.5 


[76] 


BSk4 


0.157 


-15.8 


28.0 


13.2 


14.8 


12.5 


26.4 





44.4 


-6.5 


-51.7 


[72] 


Skz3 


0.160 


-16.0 


32.0 


17.5 


14.5 


13.0 


35.1 





43.4 


10.3 


-75.8 


[7^ 


BSk8 


0.159 


-15.8 


28.0 


15.3 


12.7 


14.9 


30.6 





38.2 


7.5 


-61.4 


[79] 


Dutta 


0.162 


-16.0 


26.6 


12.4 


14.2 


16.5 


24.8 





42.7 





-51.0 


[80] 


Skz2 


0.160 


-16.0 


32.0 


17.5 


14.5 


16.8 


35.1 





43.4 


-10.7 


-50.9 


]m 


BSk6 


0.157 


-15.8 


28.0 


15.2 


12.8 


16.8 


30.4 





38.4 


6.2 


-58.2 


[72] 


BSk7 


0.157 


-15.8 


28.0 


15.2 


12.8 


18.0 


30.4 





38.4 


7.5 


-58.3 


[72] 


SKP 


0.163 


-16.0 


30.0 


12.4 


17.6 


19.6 


24.8 





52.7 


-26.1 


-31.9 


m 


BSk5 


0.157 


-15.8 


28.7 


13.2 


15.5 


21.4 


26.4 





46.5 


-8.0 


-43.5 


[72] 


Skzl 


0.160 


-16.0 


32.0 


17.5 


14.5 


27.7 


35.1 





43.4 


-24.7 


-26.1 
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TABLE I (Continued.) 



Model 


Po 




Esym{po) 


El 


E2 


L 


Li 


L2 


Ls 


Li 


L5 


Ref. 


SIIIs 


0.148 


-16.1 


32.0 


15.0 


17.0 


28.7 


29.9 





51.0 


-5.0 


-47.3 


[68] 


SKT6 


0.161 


-16.0 


30.0 


12.3 


17.6 


30.9 


24.7 





52.9 





-46.7 


[82] 


SKT7 


0.161 


-15.9 


29.5 


14.8 


14.7 


31.1 


29.6 





44.2 


-14.8 


-27.8 


[82] 


SKXm 


0.159 


-16.0 


31.2 


12.7 


18.5 


32.1 


25.3 





55.6 


-22.0 


-26.9 


m 


RATP 


0.160 


-16.0 


29.2 


18.4 


10.9 


32.4 


36.8 





32.6 


-20.5 


-16.5 


[83] 


SkSCW 


0.161 


-15.9 


30.0 


12.3 


17.7 


33.1 


24.6 





53.0 





-44.5 


[69] 


SKX 


0.155 


-16.1 


31.1 


12.1 


19.0 


33.2 


24.3 





56.9 


-23.7 


-24.3 


m 


MSk2 


0.157 


-15.8 


30.0 


11.6 


18.4 


33.4 


23.1 





55.3 


0.0 


-45.1 


[73] 


SKXce 


0.155 


-15.9 


30.1 


12.0 


18.2 


33.5 


23.9 





54.5 


-24.2 


-20.8 


m 


BSkl5 


0.159 


-16.0 


30.0 


15.3 


14.7 


33.6 


30.6 


0.0 


44.2 


-3.9 


-37.2 


[85] 


SKT8 


0.161 


-15.9 


29.9 


14.8 


15.1 


33.7 


29.6 





45.4 





-41.3 


[82] 


SKT9 


0.160 


-15.9 


29.8 


14.8 


15.0 


33.7 


29.5 





45.0 





-40.8 


[82] 


MSkl 


0.157 


-15.8 


30.0 


12.2 


17.8 


33.9 


24.3 





53.5 





-43.9 


[73] 


BSkl6 


0.159 


-16.1 


30.0 


15.3 


14.7 


34.9 


30.5 


0.0 


44.2 


-2.1 


-37.8 


[89] 


SkzO 


0.160 


-16.0 


32.0 


17.5 


14.5 


35.1 


35.1 





43.4 


-42.1 


-1.3 




Skyrmelp 


0.155 


-16.0 


29.4 


13.2 


16.1 


35.3 


26.4 





48.4 


-11.4 


-28.1 


m 


BSkl7 


0.159 


-16.1 


30.0 


15.3 


14.7 


36.3 


30.5 


0.0 


44.2 


-2.1 


-36.4 


[90] 


BSklO 


0.159 


-15.9 


30.0 


13.3 


16.7 


37.2 


26.6 





50.1 


-10.5 


-29.0 


[86] 


SGII 


0.158 


-15.6 


26.8 


15.5 


11.3 


37.6 


31.0 





33.9 


-16.1 


-11.2 


[87] 


BSkl2 


0.159 


-15.9 


30.0 


13.3 


16.7 


38.0 


26.5 





50.2 


-9.9 


-28.8 


[86] 


BSkll 


0.159 


-15.9 


30.0 


13.3 


16.7 


38.4 


26.5 





50.2 


-9.8 


-28.6 


[86] 


SLylO 


0.156 


-15.9 


32.0 


17.6 


14.3 


38.7 


35.3 





43.0 


15.4 


-55.0 


[88] 


BSkl3 


0.159 


-15.9 


30.0 


13.3 


16.7 


38.8 


26.5 





50.2 


-9.7 


-28.2 


[86] 


BSk9 


0.159 


-15.9 


30.0 


15.3 


14.7 


39.9 


30.6 





44.2 


10.8 


-45.7 


[79] 


KDE 


0.164 


-16.0 


32.0 


16.5 


15.4 


41.4 


33.1 





46.3 


12.3 


-50.2 


m 


BSkl4 


0.159 


-15.9 


30.0 


15.3 


14.7 


43.9 


30.5 


0.0 


44.2 


-2.0 


-28.8 


[92] 


SLy230a 


0.160 


-16.0 


32.0 


17.6 


14.4 


44.3 


35.2 





43.1 


32.0 


-66.1 


m 


KDEO 


0.161 


-16.1 


33.0 


17.2 


15.8 


45.2 


34.4 





47.4 


6.9 


-43.5 


m 


SLy8 


0.160 


-16.0 


31.4 


17.7 


13.8 


45.3 


35.3 





41.3 


13.7 


-45.1 


[88] 


SLy4 


0.160 


-16.0 


31.8 


17.6 


14.2 


45.4 


35.3 





42.5 


13.8 


-46.2 




SLyO 


0.161 


-16.0 


31.5 


17.6 


13.8 


45.4 


35.3 





41.5 


13.6 


-45.0 


[88] 


SLy3 


0.160 


-16.0 


32.1 


17.7 


14.4 


45.5 


35.3 





43.2 


13.8 


-46.8 


[88] 


SKMs 


0.160 


-15.8 


30.0 


15.6 


14.4 


45.8 


31.2 





43.3 


-19.4 


-9.3 


[60] 


SLy230b 


0.160 


-16.0 


32.0 


17.6 


14.4 


46.0 


35.3 





43.1 


13.9 


-46.4 


[124J 


SLy7 


0.158 


-15.9 


32.0 


17.7 


14.3 


47.2 


35.4 





42.8 


15.0 


-46.0 


m 


SLy6 


0.159 


-15.9 


32.0 


17.7 


14.2 


47.4 


35.4 





42.7 


14.6 


-45.3 


m 


SKb 


0.155 


-16.0 


23.9 


19.8 


4.1 


47.5 


39.6 





12.3 


-21.6 


17.3 


[87] 


SLy5 


0.160 


-16.0 


32.0 


17.6 


14.4 


48.3 


35.3 





43.1 


13.7 


-43.8 


m 


SLy2 


0.160 


-15.9 


32.3 


17.6 


14.7 


48.8 


35.2 





44.0 


13.6 


-44.1 


[88] 


SLyl 


0.160 


-16.0 


32.5 


17.6 


14.9 


48.8 


35.2 





44.7 


13.7 


-44.8 


[88] 


SKM 


0.160 


-15.8 


30.7 


15.6 


15.2 


49.3 


31.2 





45.5 


-18.0 


-9.3 


[93] 


SII 


0.148 


-16.0 


34.2 


20.1 


14.0 


50.0 


40.3 





42.0 


-19.2 


-13.1 


[66] 


Skzml 


0.160 


-16.0 


32.0 


17.5 


14.5 


54.1 


35.1 





43.4 


-47.9 


23.6 




SKT3 


0.161 


-15.9 


31.5 


12.3 


19.2 


55.3 


24.7 





57.5 





-26.8 


[82] 


SLy9 


0.151 


-15.8 


32.1 


17.8 


14.4 


55.4 


35.5 





43.2 


17.9 


-41.1 


[88] 


SKT3S 


0.160 


-16.0 


31.7 


12.3 


19.4 


55.9 


24.6 





58.2 





-26.9 


[82] 


SKTls 


0.160 


-16.0 


32.0 


12.3 


19.7 


56.1 


24.6 





59.2 





-27.7 


[82] 


SKT2 


0.161 


-15.9 


32.0 


12.3 


19.7 


56.2 


24.7 





59.0 





-27.5 


[82] 


SKTl 


0.161 


-16.0 


32.0 


12.3 


19.7 


56.2 


24.7 





59.1 





-27.5 


[82] 


MSkA 


0.153 


-16.0 


30.3 


15.1 


15.3 


57.2 


30.1 





45.9 


-16.4 


-2.4 




SkI6 


0.159 


-15.9 


29.9 


19.1 


10.8 


59.2 


38.2 





32.4 


23.1 


-34.4 


[95] 


MSLO 


0.160 


-16.0 


30.0 


15.4 


14.6 


60.0 


30.7 





43.9 


-13.2 


-1.5 


[19] 


SkI4 


0.160 


-15.9 


29.5 


18.9 


10.6 


60.4 


37.8 





31.7 


21.4 


-30.6 


[96] 


LNS 


0.175 


-15.3 


33.4 


15.8 


17.7 


61.5 


31.5 





53.0 


-12.8 


-10.2 


m 
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TABLE I (Continued.) 



Model 


Po 




Esym{po) 


El 


E2 


L 


Li 


L2 


Li 


L4 


is 


Ref. 


SIV 


0.151 


-16.0 


31.2 


25.1 


6.1 


63.5 


50.2 





18.4 


-23.5 


18.5 


[78] 


SGI 


0.154 


-15.9 


28.3 


19.7 


8.6 


63.9 


39.5 





25.8 


-7.0 


5.6 


m 


SKOs 


0.160 


-15.8 


31.9 


13.7 


18.2 


68.9 


27.4 





54.7 


-2.4 


-10.8 


[98] 


SkMP 


0.157 


-15.6 


29.9 


18.6 


11.3 


70.3 


37.1 





34.0 


-13.1 


12.3 


[99] 


Ska 


0.155 


-16.0 


32.9 


19.8 


13.1 


74.6 


39.6 





39.4 


-21.6 


17.3 


[100] 


SKO 


0.160 


-15.8 


32.0 


13.7 


18.2 


79.1 


27.5 





54.7 


-4.3 


1.3 


[98] 


SKYT 


0.148 


-15.4 


33.7 


19.3 


14.3 


80.8 


38.7 





43.0 


-10.8 


9.9 


\m 


Rsigma-fit 


0.158 


-15.6 


30.6 


15.5 


15.0 


85.7 


31.1 


0.0 


45.1 


-14.6 


24.1 


[67] 


SK272 


0.155 


-16.3 


37.4 


15.6 


21.8 


91.7 


31.2 





65.4 


-17.6 


12.6 


[102] 


Gsigma-fit 


0.158 


-15.6 


31.4 


15.5 


15.9 


94.0 


31.0 


0.0 


47.6 


-14.6 


30.0 


[67] 


SKT4 


0.159 


-16.0 


35.5 


12.2 


23.2 


94.1 


24.5 





69.7 








[82] 


SK255 


0.157 


-16.3 


37.4 


15.2 


22.2 


95.1 


30.5 





66.5 


-20.9 


19.0 


[102] 


SV 


0.155 


-16.0 


32.8 


31.4 


1.4 


96.1 


62.8 





4.2 


-29.1 


58.2 


[67] 


SKT5 


0.164 


-16.0 


37.0 


12.5 


24.5 


98.5 


25.0 





73.6 








[82] 


SkI3 


0.158 


-16.0 


34.8 


21.1 


13.8 


100.5 


42.1 





41.3 


35.5 


-18.4 


[96] 


SkI2 


0.158 


-15.8 


33.4 


17.7 


15.6 


104.3 


35.5 





46.9 


15.7 


6.2 


[96] 


SkI5 


0.156 


-15.8 


36.6 


20.8 


15.8 


129.3 


41.7 





47.4 


35.0 


5.2 


[96] 


Skll 


0.160 


-16.0 


37.5 


17.7 


19.8 


161.1 


35.5 





59.3 


14.2 


52.0 


[96] 



ing the empirical nucleon-nucleon scattering amp litude 
determined by Murdock and Horowitz (MH) [lQ5l | with 
isospin-dependent and isospin-independent Pauh block- 
ing corre ction s as weh as by McNeil, Ray, and Wallace 
(MRW) |l06| . the relati visti c Dirac-Brueckner-Hartrec- 
Fock (DBHF) theory [lOl\ . and the non- relativistic 
Brueckner-Hartree-Fock (BHF) theory with and without 
the t hree-body force (TBF) rearrangement contribution 
[l08l |. For these microscopic results, one can see that 
they are all consistent with each other around and be- 
low pq although there still exist larger uncertainties at 
higher density of p = 2pQ. It is interesting to see that 
the momentum dependence of the Usym,i{p, k) from the 
MDI interaction with x = are in good agreement with 
the results from the microscopic approaches. It should be 
noted that the momentum dependence of the Usym,i{Pi k) 
at po is the same for x = —1, 0, and 1 since Ugym.iip, k) 
is independent of the x parameter at po by construction 
as mentioned previously. 

Similarly as in Fig. [51 we plot in Fig. [HI and Fig. 
[TU]the the momentum dependence of the Usym,iip, k) at 
p = 0.5/90, Po and 2po using the Skyrme Hartree-Fock 
approach with MSLO, SLy4 and SKM*, and the Gogny 
Hartree-Fock approach with Dl, DIS, DIN and DIM, 
respectively. On the one hand, it is seen from Fig. [S] 
that, for < p < 700 MeV/c, the results from the MSLO 
and SKM* interactions agree well with those from mi- 
croscopic approaches while the results from the SLy4 in- 
teraction seem to display large deviation from the mi- 
croscopic results, especially at higher nucleon momenta. 
On the other hand, one can see from Fig. [10] that at 
p ~ O.bpo and po, the results from Dl interaction are con- 
sistent with those of the microscopic calculations while 
the results from DIS, DIN and DIM exhibit large dcvia- 
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0.0 0.5 LO 0.0 0.5 1.0 0.0 0.5 1.0 1.5 
p (GeV/c) 



FIG. 8: (Color online) Momentum dependence of the 
Usy,-n.i{p,k) at p — 0.5po (a), po (b) and 2po (c) using the 
MDI interaction with x = —1, 0, and 1. The corresponding 
results from several microscopic approaches are also included 
for comparison (See the text for details). 



tion from the the microscopic calculations except at low 
nucleon momenta (p < 300 MeV/c). At p = 2po, the 
Gogny Hartree-Fock calculations display different results 
from the microscopic ones and the strong model depen- 
dence appears. 

Overall, one can see from Fig. [51 Fig. [HI and Fig. 
1101 that the momentum dependence of the Usym,i varies 
from one interaction to another for the MDI interac- 
tion model, the Skyrme Hartrce-Fock approach and the 
Gogny Hartree-Fock approach. For the MDI interaction 
with X = —1, 0, and 1, the Gogny interaction with Dl, 
and the Skyrme interaction with MSLO and SKM*, the 
value of Usym,i can be negative at higher momentum 
while for the Gogny interaction with DIS, DIN, DIM and 
the Skyrme interaction with SLy4, the value of Usym,i is 
positive at higher momentum. 
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FIG. 9: (Color online) Same as Fig. [8] but for the Skyrme 
Hartree-Fock approach with MSLO, SLy4, and SKM*. 



about Usym,2 will be very useful and important for con- 
straining the theoretical models. 
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FIG. 10: (Color online) Same as Fig. [8]but for Cogny Hartree- 
Fock approach with Dl, DIS, DIN, and DIM. 



D. The second-order symmetry potential 

Usym,2iP, k) 

As for Usym,2, to our best knowledge, there is not any 
experimental/empirical information or theoretical pre- 
dictions so far. Fig. [Til and Fig. [T2l show the momentum 
dependence of Usym,2 at Q.bpo, po and 2pQ in the MDI in- 
teraction model and the Gogny Hartree-Fock approach, 
respectively. Since Usym,2 is independent of the nucleon 
momentum in the Skyrme Hartree-Fock approach, we 
only show here its density dependence in Fig. 1131 From 
Fig. [TT] and Fig.fl^ it is interesting to see that for all in- 
teractions in the MDI interaction model and the Gogny 
Hartree-Fock approach at p = O.Spoj Po and 2po, Usym,2 
firstly decreases with nucleon momentum and then es- 
sentially saturates when the nucleon momentum is larger 
than about 500 MeV/c. Especially, the results from 
the MDI interaction with x ~ 1 seem to be in quantita- 
tive agreement with those from the Gogny Hartree-Fock 
approach. Furthermore, one can see from Fig. [TT] and 
Fig. [T2| that the magnitude of Usym.2 increases with the 
density, and this is also true for the Skyrme Hartree- 
Fock approach as shown in Fig. [T31 Another interest- 
ing feature is that for the MDI interaction model and 
the Skyrme Hartree-Fock approach, Usym,2 can be ei- 
ther negative or positive while it is always negative for 
the Gogny Hartree-Fock approach with the interactions 
considered here. Therefore, any experimental constraints 



FIG. 11: (Color online) Momentum dependence of the 
Usym.2{p,k) at p = 0.5po (a), po (b) and 2po (c) using the 
MDI interaction with x = —1, 0, and 1. 
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FIG. 12: (Color online) Same as Fig. [U] but for Gogny 
Hartree-Fock approach with Dl, DIS, DIN, and DIM. 
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FIG. 13: (Color online) Density dependence of Usym,2 for the 
Skyrme Hartree-Fock approach with MSLO, SLy4, and SKM*. 

Comparing Fig.[8|with Fig.[TT]for the MDI interaction 
model, Fig. [9| with Fig. [T3| for the Skyrme Hartree-Fock 
approach, and Fig.[TU|with Fig.[T2]for the Gogny Hartree- 
Fock approach, one can see that, at fixed density and 
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momentum, the magnitude of Usym.2 is generally compa- 
rable with that of Usym,i- So, in Eq. ([S]), compared with 
Usym,i{p, k)S, the contribution from Usym,2{p, k)6'^ could 
be negligible only if S is small {S <^ 1). Therefore, we 
conclude from the present model calculations that, when 
6 is small ((5 <C 1), the Lane potential might be a good ap- 
proximation to the single- nucleon potential C/„/j,(p, S,p). 
However, the contributions of Usym,2{pTk) might not be 
simply neglected when 5 is close to 1. 



IV. SUMMARY 

Using the Hugenholtz-Van Hove theorem, we have 
explicitly and analytically expressed the symmetry en- 
ergy Esym{p) and its density slope L(p) in terms of the 
single-nucleon potential in asymmetric nuclear matter 
that might be extracted from experiments. We have care- 
fully checked the contributions of each decomposed term, 
i.e., Ei{p) and E2[p) for E^ymip), and Li(yo), i2(p), 
L^{p), L4{p) and Lq{p) for L{p) by using three popular 
phenomenological nuclear interaction models in nuclear 
structure and reaction studies, namely, the isospin and 
momentum dependent MDI model, the Skyrme Hartree- 
Fock approach, and the Gogny Hartree-Fock approach. 

Our results have indicated that the E2{p) due to the 
symmetry potential Usym,i{p, kp) is comparable with the 
Ei{p) which describes the kinetic part including the nu- 
cleon effective mass contribution and the observed dif- 
ferent density behaviors of Egym (p) for different interac- 
tions are essentially due to the variation of the symme- 
try potential Usym,i{Pi k). For the density slope parame- 
ter interestingly, we have found that, although the 
term L2 due to the momentum dependence of the nucleon 
effective mass might not have significant contributions, 
the term L5, which is from the second-order symmetry 
potential Usym.2ip,kp), generally cannot be negligible. 

By analyzing the density and momentum dependence 
of Usym,iip, k) and Usym,2{p, k) for the three nuclear ef- 
fective interaction models, we have demonstrated that 
the magnitude of the second-order symmetry potential 
Usym,2{p,k) is generally comparable with that of the 
symmetry potential Usym,i{p,k) and thus the Lane po- 
tential U„/p{p, S, k) K. Uo{p, k) ± Usym.iip, k)6 could be a 
good approximation to the single-nucleon potential only 
if the isospin asymmetry 5 is small {5 ^ 1). How- 
ever, Usym,2ip,k) might not be neglected in describing 
the single-nucleon potential Un/p{p,S,k) in extremely 
neutron(proton)-rich nuclear matter, e.g., in neutron 
stars and the neutron-skin region of heavy nuclei, where 
the value of 5 could be very large (close to 1). 

While the momentum dependence of Uo{p,k) and 
Usym,i{p,k) has been extensively investigated and rel- 
atively well constrained from the measured nucleon op- 
tical model potentials, heavy ion collisions experiments, 
and the microscopic calculations, especially around and 
below nuclear matter saturation density, knowledge on 
Usym,2{p,k) is still very poorly known. The results on 



Usym,2{p,k) presented here from the three phenomeno- 
logical models have indicated large model dependence. 
Therefore, to constrain Usym,2{p, k) from experiments or 
microscopic calculations (e.g., BHF and DBHF) based 
on nucleon-nucleon interactions derived from scattering 
phase shifts is crucial for a complete and more precise de- 
scription for Un/p{p, S, k), and thus for Esym{p) and L{p). 
Experimentally, information on the momentum depen- 
dence of Usym,2{po, k) can be in principle obtained from 
the isospin dependent nucleon optical model potentials. 
On the other hand, all analyses in the present work are 
based on the non-relativistic models, it will be thus in- 
teresting to see how our results change in relativistic co- 
variant energy-density functional, such as the relativistic 
mean field models. These studies are in progress. 
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Appendix: Models for nuclear effective interactions 

For completeness, we briefly introduce in this Ap- 
pendix the nuclear interaction models used in this work 
and also present some important expressions. These 
models include the isospin- and momentum-dependent 
MDI interaction, the Hartree-Fock approach based on 
Skyrme interactions, and the Hartree-Fock approach 
based on the finite-range Gogny interactions. These 
models have been extensively used in nuclear structure 
studies and transport model simulations for heavy ion 
collisions. 



1. Isospin- and momentum-dependent MDI 
interaction 

The isospin- and momentum-dependent MDI interac- 
tion is a phenomenological effective interacti on b ased on 
a modified finite range Gogny interaction (36lll09| . In the 
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MDI interaction, the potential energy density Spot (p, 5) 
of an asymmetric nuclear matter at total density p and 



isospin asymmetry 5 is given by 



£pot{p,5) = 



Au{x)pnPj, ^ Ai{x) 



a + 1 p^ 



Po 



dpdp' 



fAp)fr'ip') 

l + (p-p')VA2' 



(A.l) 



r 



where Au{x) 



-95.98- 



2B 



(MeV), B 
-11.70 



(MeV), Ai{x) ^ -120.57+ 
106.35 (MeV), a = 4/3, C^^r = 
(MeV), Cr,-r = -103.40 (MeV), and A = 
fi,(37r^po/2)^^'^ are obtained from fitting the momen- 
tum dependence of single-nucleon potential to that pre- 
dicted by the Gogny Hartree-Fock and/or the Brueckner- 
Hartree-Fock calculations, the saturation properties of 
symmetric nuclear matter and the symmetry energy of 
30.5 MeV at nuclear matter saturation density pq = 0.16 
fm~'^. The incompressibility for cold symmetric nuclear 
matter at saturation density po is set to be Kq = 211 
MeV. The x parameter in the MDI interaction is intro- 
duced to vary the density dependence of the nuclear sym- 
metry energy while keepin g ot her properties of the nu- 
clear equation of state fixed |36| , and it can be adjusted to 
mimic the predictions of microscopic and/or phenomeno- 



logical many-body theories on the density dependence 
of nuclear matter symmetry energy. We would like to 
point out that the MDI interaction has been extensively 
used in the transport model for studying isospin effects 
in intermediate energy hea.vy-io n collisions induced by 
neutron- rich nuclei [H, Ill0l4ll7j | , the study o f the ther - 
mal properties of asymmetric nucl ear matte r jllSl Ill9l | , 
and the compact star physics [l^, Il20l . Il2l| . In partic- 
ular, the isospin diffusion data from NSCL/MSU have 
constrained the value of x to between and —1 for nu- 
clear matter densities less than about 1.2po [36j . 

In mean-field approximation, the single-nucleon po- 
tential for a nucleon with momentum p and isospin 
T in asy mmetric nuclear matter can be expressed as 
p9l[T22l . lT23j 



J 



Ur{p,6,p) ^ Au{x)^ + Ai{x)^ 
Po Po 



B 



Po a + 1 p^ 



r-l 



2C7- 2 

Po 

2 C7- . — -7- 2 



PO 



■ttA' 



2M . 

In- r- 



Po 

— : — /n ^^ H — — 2arctan ^ PFr _^ 2arctan- — ^ 

{p-pF^y+A-^ AAA 

A^ 2pF_^ P + PF-t P ^ PF-^ 

' ^ 2arctan : + 2arctan : 

A A 



where p^ = p(l + tS)/2 and pF^ ~ h(3Tr'^pr) 



A '' "''7 

2pA (p — pf. 

,1/3 



■_J2 + A2 



A 



,(A.2) 



2. Skyrme Hartree-Fock approach 



For the Skyrme interaction, we use the standard form jl24| that has been shown to be very successful in describ- 
ing the structure of finite nuclei. Neglecting the spin-orbit interaction term which is irrelevant in nuclear matter 
calculations considered here, the nuclear effective int eract ion in the standard Skyrme interaction is taken to have a 
zero-range, density- and momentum-dependent form |124| . i.e.. 



Skyrme 



1 
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(ri, ra) = t^{l + xoPa)5{r) + -h{l + xiP„)[P'^5{r) + 5{r)P''-] +t2{l + X2Pa)P' ■ 5{r)P 



-H{l + x^PMR)rKr), 
6 



(A.3) 



where r = r2, R = (ri + r2)/2, P„ is spin exchange operator, P = ^(Vi — V2) is relative momentum operator 
acting on the right and P' is its conjugate which acts on the left. The to, xq, ti, Xi, t2, X2, ^3, X3 and a are the 9 
Skyrme interaction parameters determined from fitting the binding energies, charge radii, and other properties of a 
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large number of nuclei in the periodic table. In the Skyrmc Hartree-Fock approach, the single-nucleon potential in 
asymmetric nuclear matter is given by jl24l | 

P P 1 

Ur{p,5, k) = —p[ti{2 + xi) + <2(2 + 2-2)] + —PrMl + 22:2) -h{l + 2xi)] + -to[(2 + x^)p - {2X0 + l)Pr] 

+ ^i3p"[(2 + X3)p - (2x-3 + l)pr] + ^i3p"-M(2 + X3)p' - {2x3 + l){pl + P^)] 

+ i[ii(2 + X,) + t,{2 + -2)](^ + ^) + l[h{2x, + 1) - t,{2x, + 1)]^. (A.4) 



3. Gogny interaction 



Gogn y intera ction has been proved to be very successful in describing not only nuclear structure but also nuclear 
matter |62l . Il25j . Compared with the Skyrme interaction which only contains 6 forces, Gogny interaction is featured 
by two finite range terms plus one S force that can well mimic the nucleon-nucleo n effe ctive interaction. Neglecting 
the spin-orbit interaction term, the conventional Gogny interaction is given by (62l . Il25j l 



T/if ^"^(n, r2) = ^(W, + B,P, - H,P^ - M,P,Pr)e~^"^f^' + to{l + xoP,)[p{R)rS{r), 



(A.5) 



where Wi, Bi, Hi, A/i, pi, W2, B2, H2, A/2, p2, ^o, 2:0 and a are the 13 Gogny interaction parameters, and Pt is the 
isospin exchange operator. By using Hartree-Fock approach, we explicitly write down its single-nucleon potential 



2 2 
Ur{p,S,k) ^ pY^^%^^^{w, + ^)-p^Y.^i^cl{H, + ^) + top^^ 



i=l 



(1 + y)^" {\+Xo)Pr 



1 

-toap"+M3 - (2x0 + + y Z,{k, t){-W, - 2B, + H, + 2Af,) 
1 

yZ^{k,-T){H, + 2M^), 



(A.6) 



with 



Zi{k,T) 



y/npik 

where erf{x) is error function. 



+ l{erf[p^{k + fci.J/2] - erf[p,{k - fcFj/2]}, 



(A.7) 
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